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Q-r Introduction 

It is well known [1] that investigation of the problem occurrence of chaotic 
dynamics naturally led to the need to study simple models in which there is 
chaos. Further study of the real physical processes in which there is chaos 
showed that many simple models reflect its main, for example, the occurrence 
>- ' of stochastic separatrix layer, formation of nonlinear resonances and etc. 

One-dimensional chains of particles have also been the subject of research 
■^J- ■ in studying the problems of chaotic dynamics [1] because some problems for 

them are reduced to the study of standard maps known for nonlinear dynamic 
systems. On the other hand investigation of the phenomena of chaos chain 
of particles is also important in connection with the active use of the particle 
method for simulating the behavior of materials. Speaking about this method 
we consider the particles as point masses, and not as discrete elements allow- 
ing to reduce the equations of continuum mechanics to the difference system 
of ordinary differential equations. One of the most well-developed variants 
of the particle method is the method of molecular dynamics in the classical 
version of which particles act as atoms and molecules. If interatomic poten- 
tial is known then the dynamics of molecular compounds can be modeled 
with high accuracy. 

Despite the good correlation between the results of computer simulation 
observed experiment and material behavior there is a problem of under- 
standing their internal mechanisms in terms of nonlinear dynamics. For one- 
dimensional system corresponding model problem is reduced to consideration 
of the particle mass m in the cell. The dynamics of particle is determined by 
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the Hamiltonian 
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H = H(p, q, t) = |- + V(q), V(q) = U(q) + U(£(t) - q), (1) 

in which p is the particle momentum, and form of the function £ (t) depends 
on the collective behavior of particles in the system. Formula (PQ) means that 
particle interacts with a stationary particle on the left side and interacts 
with a particle moving in a law q2 = £(£) on the right one. The interaction 
between the particles is characterized by the potential U(q). 

The phenomenon of strong chaos for particles interacting by means of 
the Lennard- Jones potential U(q) was considered in [2]. Chirikov criterion of 
resonance overlapping [3] was used as the condition of chaos. The function 
K Hjn -i corresponding to the criterion was calculated numerically with respect 
to energy. In this article we obtain the asymptotic representation for this 
function at low and high energies. 

1 Chirikov criterion of resonance overlapping 

Let us remind the model [2] . Interaction is determined by the Lennard- Jones 



potential U(q) = U 



a \ a l a 



with a > j3 > 0. The perturbation 
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£(t) is given by the formula £(t) = 2a (l + e(t)), e(t) = aicoswit with 
| ckx j 1. We introduce the dimensionless variables setting H — > H/Uq, 
q — > q/a — 1, t — >■ t/to, 0J\ — > oJit^ where t = ^Jma 2 /U$. Then Hamiltonian 
([1]) is written in the form 
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H(p, Q,t) = J + V(q), V(q) = U{q + 1) + U{1 + e(t) - q). (2) 

Given the small amplitude e(t) we carry out the expansion of the Hamiltonian 
(J2J) with respect to e(t): 



H(p,q,t) = H (p,q) - e{t)H x {q) + 



V 
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H (p,q) = ^+W(q), H x {q) = ^(1-g), W(q) = U(q+1)+U(q-1). (3) 

Action-angle variables (/, ip) are commonly used in studying the dynamics 
of Hamiltonian systems [1]. The unperturbed action is conserved along the 
phase trajectory and its value is determined by initial conditions (po,qo). If 
you choose po = the correspondence between the unperturbed action / and 
the parameter go is unique. Transformation to the action-angle variables can 
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be carried out numerically in accordance with formulae 
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HQo) = ^ J p(qo,q)dq, <p(q,q ) = w J p ^ q y u = qT~ ' ( 4 ) 



go 



p(qo, q) 



where p(qo,q) = y 2(V(g ) — V(q)). We express the variables (p,q) in the 
terms of Then Hamiltonian H(p,q,t) can be expanded in a Fourier 

series: 



H(p,q,t) = ff (I)-^«i y^if n (/)[cos(wy+y -a;it)+cos(wy+yo+a;it)], 

n 

Hn = lj d V U^\l-q) 
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n 

2n (5) 

cos nip. 



Let us consider the resonance condition for a three of numbers (n,I n ,ui): 
nu{I n ) = wi. 

To analyze the dynamics of the system in the vicinity of nonlinear reso- 
nance we use the standard method [1]. In this case the only resonant har- 
monic in fl5]) is fixed: 

H(p, q, t) = H (I) - -a\H n (I) cos(ny2+y? -wit). (6) 

We introduce a new phase ip = rup + ip — uj\t then the equations of motion 
corresponding to Hamiltonian fl^]) have the form 

dl dH 1 . 

at 2 

dip dH 1 d 

— = — = «,(/) - a* - -ax cos^-tf n (/). 

As the equations (J7J) are examined in the vicinity of the resonance action I n \ 
it is assumed that the value J = I—I n is small. We expand H (I), and 
neglect the terms J 2 correspondingly. Then the equations are written in 
the form: 

The obtained equations correspond to the so-called nonlinear pendulum Hamil- 
tonian [1]: 

H(J,ip) = ~~ -aiH n (I n ) cosi/j. 
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The width of the resonance AJis calculated from the condition: 



In terms of frequency u we have 



Au^^-a.HM^iQ. (8) 

Consider the invariant curves for the resonant action I n , I n -i- In order to 
form a chaotic region in the plane (p, q) it is necessary to destroy the curves 
corresponding to these values of the action. Let us assume that AQ (AI) 
is the corresponding width of the resonance and 5Q = co(I n ) — u(Ir n -i)) 
(8I=I n — I(n-i)) is the distance between resonances. The parameter K char- 
acterizing the degree of resonance overlapping is equal to K=AI / 5I~ Af2 / 8£l 
[1]. In accordance to Chirikov criterion [3] the overlap of resonance takes 
place on condition that \K\ > 1 : 

V ( ( \\ A ^( J n) + AU (/(„_!)) 

^n,«-l(5b(Wi)) = ry r 1. (9) 

In [2] the function K 32 was calculated numerically at different \a%\ <C 1 
with respect to q (energy). In this article we construct an asymptotic rep- 
resentation for the functions Ki^(qo(coi)) , K 2 ,z (go(^i)) at small and large 
energies. 



2 Approximate action-angle variables at small energy 

The criterion (j5j) includes the frequency and its derivative (|SJ), and the value 
of the Fourier coefficient (j^J). From the point of view of physics the smallness 
of the energy values correspond to the fulfillment of the inequality £< 1, 
This allows us to use an approximation for W(q) ([2D as polynomial functions 
with respect to q and construct approximate formulae for the action-angle 
variables at E — > 0. Approximation for the potential W is obtained after 
expansion in the vicinity of the point q = 0: 

W app {q) = W + bq 2 + cq\ b = U®(l), c = (10) 
Let us introduce the approximate action-angle variables assuming 



W2 q 
I f f 

Qi Qi 



dq 

Papp 
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where p app = y2(E — W app (q)), and the coordinates of points Qi are deter- 

d E 

mined from the equation E = W app {Qi). The quantity Q = — — = 2 — is 

oJ TyE) 

the cyclic frequency where T(E) is the period of motion in the potential (fTUl) . 
Omitting simple calculations we have the following presentation for T(E): 

Q2 Ti 

T f dq If dl 
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[ ^2(E-W app (q)) V2-c_J Ti y/W-PXAt + P) 



Substituting t = B cos ip we write this expression in the form : 

tt/2 

T 12 1 f dt 12 1 

K(k), 11 



2 V c ^Tf + A 2 7 Vl-^sinV V c ^Tf + A 2 
where K(fc) is the complete elliptic integral. As a result obtain 

where fc) is the elliptic integral of the first order. 

3 Asymptotic behavior of different quantities at small 
energy 

dT dT 

We are interested in the asymptotic behavior of T and the derivative — -=— — 

on oE 

at H — y 0, then 

* = ^(i + £*<*)), -^L=^(i-f) + o(^). (13) 

The smallness of allows us to use the asymptotic formula for K(fc): 

K «-l{ 1+ @ 2 * 2 + (^) 2fc4 + -^ < 14 > 

From here and (fTTT) we have 

l=*(l-**£)+0{H*). 
2 ^/2b\ ) ~ y > 
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Hence we immediately obtain the leading order in the expression for the 
frequency and its derivative which can be written as: 



2b + 0(H 2 



2UW(1) + 0(H 2 ) 



dtu _ dtu 3c n /TT\ _ 



+ 0(H). 



(15) 



Let us write the function q(J,ip) in the terms of variables (J, ip) . The 
error does not exceed 0(H). Because of k ~ \/H (TT5|) we use asymptotic 
formulae for K(k) ([SD and k) = $ + ... It results in ip = § + 0(H) and 
$ = if + 0(H). 

To calculate the Fourier coefficients (J5J) we expand the potential in the 
vicinity of q = 0: 



2tt 



7T 



i/ n = — dtp cos mp 



uM{l) q +U® (1)1 + 



The function g( J, r/>) is equal to q — T\ cos((p+0(H)) where T\ = a/ H/b + 
Q_(H 2 ). Then we have 



H 9 ~ H 



U®(1) 
4f/(2)(l)' 



From dSJ), ([I5D we have inequalities |Aw(ii)| ~ \JH > |Aw(J 2 )| ~ vH > 
|Ao;(/3)| ~ \J H 3 . It allows us to obtain expression for functions | i^2,i |, |-^3,2| : 




Aw(/i 



Aw(I 2 



cut y 2 V v ; 8t/(2)(i) 



^ Wl \/2 ai 32£;( 2 ) (1) f/(2)(l)- 



(16) 



4 Approximate action-angle variables at high energy 
and asymptotic behavior of different quantities 

From the point of view of physics high values of energy E correspond to the 
fulfillment of the inequality E/Uq^>1. Since the potential W(q) = U(l+q) + 
U(l—q) is symmetric function with respect to q — > —q it is enough to consider 
it at the interval (0, 1). It is clear that for large values of E turning point 
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coordinate qi ~ 1. Then the leading contribution in V(q) is determined by 
£7(1 — q). We introduce an approximate potential 



WuM 



:i-<?) c 



and action- angle variables Ju g , > 



(17) 



Jbig{E) 
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Pbigdq, 



^biq^b 



big ^ big j 



dE 



big 



dJ, 



big 



'i<J 



Qi 



dq 

Pbig'' 
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where function p big = y2(E — Wu g (q)), and the coordinate Qi is defined by 
the equation 

E = W big (Q 1 ). (19) 

From a mathematical point of view, values /, Ju g determine the area in phase 
space for functions p, pbi g respectively. Area difference can be estimated by 
the sum of functions 

i 

\I(E) - J b ig(E)\ < ~[\Qi - qi\ max|p| + | I \p - Pbig ]dq\). (20) 

Qi 

Leading order of Q\— q± is determined from ( TTTj) . (lT9~j) and is equal to Q\—q\ = 

— 3— — (l + 0(Qi ?_/3 ))- Since max \p\ < y/2E then the first term in (T2"0"j) has 

y 

3 ft 

is equal to 



the following order (Qi — q\)y E ~ . . The integral term in f[2"0"j) 

E L/A+(L Pj/Oi 



1 








Qi 





Whig ~ W 
P + Ph 



dq 



Qi 



it-) 



Hence 



1(E) - J big {E) 



dq 

Pbig 



1 



Qi 



l-p+a/2 



1 



E l/2+(X-p)/a- 



(21) 



E X/2+(l-0)/a 

on condition that a > 2(/3 — 1) which provides a decrease in the right side 
of f)2ip . Differentiation fl2T]) with respect to results in relations for the 
frequencies corresponding to the action and their derivatives: 



du dfl 



big 



l-p)/a ' Q E Q E 
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E 5/2+(l-0)/a ■ 



(22) 



The difference between the phases cp and is written in the form: 

\(f — 0| < \bj — Vt bi g\tbig + U)\t — thig\ ~ £3/2+(l-f3)/a U \ U ~~ ^ bi 9 1 



For high energy frequency u ~ ~ v -E then from here and (1221) we have 

Let us transform from phase to phase in the integral (EJ with accuracy 
( 123]) and use the equality q((p + it, Jbi g ) = Qi — q((fi, Jug)- Since 2 = Q 2 + Qi 
we rewrite fl5]) in the following form: 

H n = \ fdq [uW(2-q) + U^(Q 1 + q )(-ir] C0S -^^-q. 
Qi 

A leading term with respect to E is determined by means of integration: 



I 7VT> 

H n = - [-U(L -q) + U(Q 1 + q)(-ir) cos — q 

Q> 



Qi 



1 7T f TTT) 

2 0^01 J lUiL - q) - U(Ql +q){ - ir] 8i " Q^Q?- (24) 



The outside terms have a singular behavior with respect to E since 
(-1)" , , „, „ v , f-l) n 1 



-17(30 - U(2Q 1 )} = - ^^[1 + 1/21 + O (Vj 



The integral in fl2^|) is smaller than these terms then the Fourier coefficient 
is equal to 

#„=-#fci2[l + l/2«] + ... (25) 
From Q, fl2"5"j) we obtain expression for functions K nn _i : 



Kn.n-1 = y y(l + l/2 Q )(n - 1/2). 

The parameter a = 12 for the Lennard- Jones potential then we obtain 
K 2A = 1.5^/2, ^3,2 = 2.5v^I72. 
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5 Comparison with numerical results for ^2,1(90) 

The function K 2 ,i(9o) was defined numerically (exact) and calculated (asympt) 
in accordance with formula (1161) . Its behavior depending on the initial posi- 
tion of the resonance trajectory is presented in Fig.l (exact). From here it 
is seen that formula (TTB]) is a good approximation for K 2 ,i(qo) as for small as 
high energies. 

K 2 ,i , 




1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

0.05 0.1 0.15 0.2 q 

Fig. 1. Graph of the function ^2,1 (go) defined numerically (exact) and 
calculated (asympt) in accordance with formula (Tl6l) . 
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